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Abstract 

We consider Einstein's equations coupled to the Euler equations 
in plane symmetry, with compact spatial slices and constant mean 
curvature time. We show that for a wide variety of equations of state 
and a large class of initial data, classical solutions break down in 
finite time. The key mathematical result is a new theorem on the 
breakdown of solutions of systems of balance laws. We also show that 
an extension of the solution is possible if the spatial derivatives of 
the energy density and the velocity are bounded, indicating that the 
breakdown is really due to the formation of shock waves. 



1 Introduction 

A question of central interest in general relativity is that of the long-time 
behaviour of self-gravitating matter. Mathematically the starting point is 
to get suitable existence and uniqueness theorems for the Einstein equations 
coupled to the equations of motion of the matter. Since the resulting system 
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of partial differential equations is difficult to handle it makes sense to begin 
with solutions of high symmetry. One of the most popular matter models in 
applications is the perfect fluid described by the relativistic Euler equations. 
It is to be expected that, as in classical hydrodynamics, a major difficulty in 
studying the long-time behaviour of solutions of these equations is the for- 
mation of shocks from smooth initial data. If the solution is to be continued 
beyond these it is necessary to leave the realm of classical solutions of the 
equations. It should be noted that until very recently most of the theorems 
on solutions of the Euler equations involving shocks were a one-dimensional 
context (plane symmetric solutions). This may change following the recent 
work of Christodoulou on the formation of shocks in special relativity without 
symmetry assumptions. 

One way to avoid the difficulties involved with fluids is to consider instead 
collisionless matter described by the Vlasov equation. In that case quite 
a lot of mathematical results are available for the coupled Einstein-matter 
equations [15, 2]. If, on the other hand, we face the problems associated 
with the fluid description, as we do in this paper, several natural questions 
arise. First, are there global existence theorems for weak solutions of the 
special relativistic Euler equations? Some positive answers have been given 
in [17], [4] and [7]. Second, do these results extend to the case of a self- 
gravitating fluid? Theorems have been proved under certain assumptions by 
[8] and by [3]. Third, can it be proved that classical solutions of the Einstein- 
Euler equations break down in finite time? It is the third question which is 
addressed in the following. 

The symmetry assumed in this paper is plane symmetry where the solu- 
tions are invariant under the action of the full isometry group of the Euclidean 
plane. In particular this reduces the full problem to a problem in one time 
and one space dimension. It is assumed that the position space is compact. 
This circumvents the need for boundary conditions and is the analogue of 
periodic boundary conditions in the non-relativistic case. The notion of 'fi- 
nite time' breakdown is subtle in general relativity. In that theory there is a 
free choice of time variable and if it has been proved that a solution breaks 
down after a finite amount of a particular time coordinate it is necessary to 
think carefully about what this means geometrically. This is what is rele- 
vant for physics. In this paper we use a constant mean curvature (CMC) 
time coordinate. This means that the value of t at some point of spacetime 
is equal to the mean curvature of the unique compact CMC hypersurface 
passing through that point. 
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The main result of this paper (Theorem 5.2) shows that plane symmetric 
classical solutions of the Einstein-Euler equations exhibit finite-time break- 
down for a wide variety of equations of state. The key new analytical result 
used to prove this is a theorem on the breakdown of solutions of systems 
of balance laws which is of interest in its own right. This is combined with 
general estimates for the Einstein equations to give the final result. 

In the last section we put the result into context, considering several 
issues. The Euler equations are compared with other matter models. The 
assumptions on the equation of state are compared with cases known in 
the literature, and the relations to existing results on weak solutions are 
discussed. 

2 The Einstein equations 

Let (M, g) be a spacetime, where the manifold is assumed to be M = I x T 3 , 
I is a real interval and T 3 = S 1 x S 1 x S 1 is the three-torus. We require 
that the metric g and the matter fields are invariant under the action of 
the Euclidean group E 2 on the universal cover, and that the spacetime has 
an ^-invariant Cauchy surface of constant mean curvature (CMC). As was 
shown in [12], there is a local in time 3 + 1 decomposition of (M,g) where 
each spatial slice has constant mean curvature. We can introduce spatial 
coordinates x, y and z on each slice, with ranges [0, 2ir] and period 2-7T, and 
a CMC time coordinate t — tr k where k is the second fundamental form of 
the slices. 

The metric can be expressed as 



where a, A and /3 depend on t and x and a depends on t only. Here a is 
the lapse and (35^1 is the shift. The coordinates can be chosen such that 



ds 2 = -a 2 dt 2 + A 2 [(dx + pdt) 2 + a 2 {dy 2 + d^ 2 )] , 



(1) 



J 2n (3dx = 0. 



It is convenient to introduce the orthonormal frame 



e = a d t — a [3d, 
ei = A^d x , 
e 2 = (aA)~ 1 d y , 
e 3 = (aA)" 1 ^. 



•x ■ 



(2) 
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We will use indices a, b, . . . to denote spatial coordinate components and 
i,j, ... to denote spatial frame components. 

The frame components of the second fundamental form may be written 

kij = -\{K - t)Sij + \{W - t)5vi5^ (3) 

where K is a function of t and x. We also introduce the notation p = T 00 , 
j = T 10 and S lk = T lk for frame components of the energy momentum tensor 
T^ u . The Einstein equations can then be written 

K' + (3K - t)A' x A' = SnAj, (4a) 

(VA)" = -\A^ 2 [K 2 + l(K-t) 2 -t 2 + Wnp] , (4b) 

a" + A~ l A' a' = aA 2 [K 2 + \{K - tf + 4tt( P + tr S)] - A 2 , (4c) 

d ta = a[-f3' + ±a(3K-t)], (4d) 

d t A = -aKA + (A[3)', (4e) 
d t K = I3K' - A~ 2 a" + A~ 3 A'a' 

Ui) 

+ a[-2A- 3 A" + 2A- 4 (A') 2 + Kt-4n(2S n -trS + p)]. V ; 

Here a prime denotes differentiation by x. Equation (4a) and (4b) are the 
momentum and Hamiltonian constraints, respectively. The lapse equation 
(4c) comes from the constant mean curvature condition, while (4d) and (4e) 
are consequences of the choice of spatial coordinate conditions. The last 
equation (4f) is the only independent Einstein evolution equation in this 
case. 

Using (4) it is possible to show that many of the fundamental quantities 
are bounded. This was done in [12] for plane, hyperbolic and spherical sym- 
metry, when the matter satisfies the dominant energy and the non-negative 
pressures conditions. The non-negative pressures condition was subsequently 
relaxed to the strong energy condition in [14] (note that plane symmetry is a 
special case of local U{\) x U{\) symmetry). Also, a bound for a" 1 in terms 
of the previously obtained geometric bounds was found in [13]. 

Theorem 2.1 (See [12], [14] and [13]). Let a solution of the Einstein equa- 
tions with plane symmetry be given and suppose that when coordinates are 
chosen which cast the metric into the form (1) with constant mean curvature 
time slices the time coordinate takes all values in the finite interval \to,t\) 
with ti < 0. Suppose further that the dominant and strong energy conditions 
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hold. Then the following quantities are bounded on the interval [to,ti): 

a, a" 1 , a', A, A' 1 , A', d t A, a, a" 1 , d t a, K, (3, (3'. (5) 
The bounds involve only 

t ,h, ||a||o, lla^Ho, Pllo, ||^ _1 ||o, (6) 
where ||-|| is the supremum norm on the initial surface t = t . 

3 The Euler equations 

In this section we will rewrite the matter equations for a perfect fluid in plane 
symmetry as a system of balance laws. 

The basic matter variables are the pressure p, the energy density /i and 
the unit 4- velocity £/ M of the fluid. Because of the plane symmetry the frame 
components of U can be written as 

U° = (l-u 2 )- 1/2 , U l = u(l - u 2 )- 1 ' 2 and U 2 = U 3 = 0, (7) 

where uG (-1,1). The energy momentum tensor for a perfect fluid is 

= {fi + p)U^U v + pg^, (8) 



which implies 



1-u 2 
■oj _ (/' • /')" 



P = T 00 = ^-p, (9a) 



11 _ mil _ ( ^+P)U 2 
1-U 2 

22 = ^33 = T 22 = T 33 = p (Qd) 



oil rpii KH' -r F " . / n x 

S =T =— —+P, (9c) 



It will be convenient to introduce two new variables 



„-(£)», (10 a) 
J \am/ m + p 
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where m is a dummy integration variable and the constant is arbitrary. 

The matter equations are given by the vanishing of the divergence ^J v T va 
of the energy momentum tensor. The spatial frame components V u T u2 and 
V u T u3 vanish identically. Expressing the remaining two components in terms 
of <f, w and u gives 



w [2ue (u) + (1 + -u 2 )ei(-u)] 

+ (1 - u 2 ) [(1 + u 2 w 2 )e (<p) + u(l + w 2 ) ei (<p)] 
- w(l - u 2 ) [Ku 2 + t - 2uA~ 1 e 1 (A) - 2ua~ 1 e 1 (a)] 



0, (11a) 



w 



[(1 + u 2 )e (u) + 2ue 1 ( 

+ (1 - u 2 ) [u(l + w 2 )e (ip) + (u 2 + w^e^ip)] 
1 - u 2 ) [{K + t)u - 2u 2 A- 1 e 1 {A) - (1 + u^a^e^a)] = 0. (lib) 



w( 



Adding and subtracting (11a) and (lib) we get 

E = w[e (u) + ei(u)] + (1 - u) [(1 + uw 2 )e ((f) + (u + u> 2 )ei(y?)] 

— w(l — u) [Ku + t- 2uA~ 1 e 1 ( y A) - (1 + u)a -1 ei(a)] = 0, (12a) 



Ei = w[e (u) - ei(u)] - (1 + u) [(1 - wu> 2 )e (y?) + (u - w 2 )ei(y?)] 

- w(l + u)[Ku-t + 2uA- 1 e 1 (A) - (1 - u)a _1 ei(a)] = 0. (12b) 

The linear combinations (1 + u)(l ± w)E + (1 — u)(l =F w)E 1 give 

'Ku + tw- 2uwA' 1 e 1 (A) 



D + u + (1 - u 2 )D + ip = (1 - u j 



D_u - (1 - u 2 )D_(p = (1 - w 2 ) 



with differentiation operators 



1 + uw 

Ku-tw + 2uwA' 1 ei(A) 
1 — uw 



— a 1 ei(a) 



(13a) 



a 1 e 1 (a) 



(13b) 



u + w u — w . . 

A- = e o + — ei and D- = e + ei. (14) 



1 + uw 



1 — -uw 
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Next, we introduce the variables 

, 1 n 1 + U , . 

r = if H — in and s = ip in . (15) 

Y 2 1-u T 2 1-u y ' 

These are analogues of the Riemann invariants found by Taub [18]. In con- 
trast with the special relativistic case, r and s are not invariant. In particular, 
it follows from (13) that instead of a system of conservation laws, r and s 
satisfy a system of balance laws: 

tw + Ku - 2uwA~ 1 e 1 (A) , . . . 

D + r = - alexia , 16a 

1 + uw 

£Ls = ^ + a _1 ei«. 16b 

1 — -uw 

Here w and w should be regarded as functions of r and s given by (15) and 
(10). 

We will now make some assumptions on the equation of state to ensure 
that the maps between the different matter variables introduced in this sec- 
tion are well behaved. 

Lemma 3.1. Suppose that the equation of state is given as a smooth function 
p of ii such that 

1. the strong energy condition (/i + p > and fi + 3p > 0) holds, 

2. the dominant energy condition (—/i < p < /i) holds, 

3. < J < 1 and 

4- 0>o- 

Then the maps 

(r,s)-( u ,<p), R 2 ^(-l,l)xR, (17) 
[/,i—>(p, (0, oo) — > (</?_, oo), — oo < </?_ < (18) 

are smooth bijections, 

(fx,u)^(p,j,S) (19) 
smooth and 1-1 on (0, oo) x (—1, 1), 

|ihmj, (0, oo) -> (0, 1) (20) 
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is smooth and increasing, and 



Proof. From (15), 



r ~\- s r — 5 

ip = — - — and u = tanh — - — , (22) 

which establishes (17) without any restrictions on the equation of state. Ap- 
plying condition 1 and 3 to (10b) shows that ip is a smooth strictly increasing 
function of /i and thus 1-1. Moreover, from condition 4 and 2 it follows that 

lim (f = oo. (23) 

Thus p is a smooth bijection from (0, oo) to (</?_, oo), where ip- = lim M _> + p> > 
— oo. 

That p, j and S are smooth functions of /i and u is evident from (9). The 
map is 1-1 since p is 1-1 because of condition 3. Finally, (20) follows directly 
from condition 3 and 4, and 

dw dw /dp 1. /dp\- 1 d 2 p , , 

^ = d^/dH<" +p >(dP (24) 

because of condition 1, 3 and 4. □ 



4 Balance laws 

In this section we will obtain a blowup result for the system of balance laws 

D + r = e (r) + re(r, s)ei(r) = f(t, x, r, s), (25a) 
D_s = e (s) + A(r, s)e 1 (s) = g(t, x, r, s), (25b) 

where {e ,ei} is a pseudo-orthonormal frame with respect to a Lorentzian 
metric on [to,ti) x S 1 for some real interval [t ? The operators D + = 
e + nei and D_ = eo + Ae x are called characteristic derivatives and we denote 
the corresponding integral curves (or characteristics) by g x and a x . The suffix 
shows where the curves intersect the surface t = t , i.e., g x (0) = (t ,x) and 
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<7x(0) = (to,x). We also denote the commutation coefficients by c\ so that 

Let be an open subset of R 2 , to be further specified below. We assume 
that at least one of k and A is genuinely nonlinear on Us, i.e., dn/dr 7^ or 
dX/ds 7^ 0. Without loss of generality we may assume that dn/dr > 0. We 
also assume that k — A 7^ on U§ so that the system is strictly hyperbolic. 

Let the initial values of r and s be r (x) = r(t ,x) and s (x) = s(to,x). 
We will show that for certain choices of tq and sq, r and s are bounded 
while d x r or d x s blow up in finite time, indicating the presence of a shock 
wave. There are similar results for systems of conservation laws (see, e.g., 
[6]), but we cannot hope for the same generality here. For more specific 
results, a better control of the source terms / and g is needed. For example, 
if the sources are superlinear in r and s we can expect blowup of r and/or s 
themselves [11]. On the other hand, global existence of smooth solutions have 
been shown under certain conditions [20]. We will leave these considerations 
aside and aim for a more general but less sharp result. 

4.1 Bounds on r and s 

First we need to specify the set Us. The initial data (r , s ) is a smooth map 
S 1 — > R 2 . The image is a plane curve with convex hull U . For any 5 > we 
define 

U S = {P e R 2 ; d(P,U ) < 5} and Q s = Kh) x S 1 x U s , (26) 

where d is the Euclidean distance in R 2 . Note that Us is convex as well, as 
follows easily from the triangle inequality. 

To be able to perform integrations along lines of constant r in a well 
defined way we need to introduce the following construction. Since S 1 is 
compact, there are points (r_,s„) and (r + , s + ) on the initial data curve in 
R 2 where r attains its minimum and maximum, respectively. These are 
also global extremal points of r in Uq since Uq is convex. By construction, 
(r_ — 5, S-) and (r + + 5, s+) are extremal points of r in the closure of Us, and 
we let s = 7(7") be the straight line between those two points. We can now 
state the following lemma. Here ||-|| is the supremum norm on S 1 . 

Lemma 4.1. The line from (r, s) to (r, 7(7")) is contained inUs for all (r, s) G 
Us- Moreover, \~f'\ < 5 _1 ||so|| and \s — 7(r)| < 2||s || + S on Us- 



9 



Proof. Since r_ — 5 < r + + S, the line s = 7(7") intersects the balls with 
radius 5 around (r_,s_) and (r + ,s + ). It follows from the convexity of Us 
that s = 7(r) is contained in except for its endpoints (r_ — 5,s-) and 
(r + + 5, s+). The first statement then follows directly from the convexity of 

The second statement is just the rough estimate 

W\ = !f+Ziz! < J!M (27) 

and the third follows from the fact that |7(V)| < ||s || and \s\ < \\s \\ +5 on 
U s . □ 

Let 

F = sup{|/|, \g\} and E = inf{e ° + ne®, e ° + Aei°}. (28) 

n s n 5 

The following lemma provides a crude estimate of r and s. 

Lemma 4.2. // F finite and E > 0, any smooth solution (r, s) of (25) 
with initial data (r , s ) remains within U s fort< min{t + 2- 1 / 2 5EF~\ tj. 

Proof. From (25) we have that |r — vq\ < Ft; along a characteristic g x as long 
as t < ti, where ^ is the parameter along g x , chosen such that D + = d/d£ 
and g x (t ) = (t ,x). By the definition of D + , dt/d£ = e ° + ne\ > E, 
so ^ < E _1 (t — t ) and \r — r \ < FE~ l {t — t ). We can obtain similar 
estimates for s along the characteristics o~ x . Thus the distance between (r, s) 
and (r , sq) is at most \/2FE~ l {t — t ), and the conclusion follows. □ 

As we saw in Lemma 3.1, [/,*—* (p: (0, oo) — > (</?_, oo) is smooth and 1-1. 
It is quite possible that </?_ is finite, i.e., that /x = for finite r and s. This is 
indeed the case for a relativistic polytropic perfect fluid, for example. This 
complication can be circumvented by a further restriction on 5. 

Lemma 4.3. Let (p = ^mi xeS i{r (x) + s (x)}. If 

< 5 < ipo - <p- (29) 

then <f > (/?_ on Us. 
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Proof. If /i = at some point on the initial surface, then tp = (r +So)/2 = P- 
there. Thus the right hand side of (29) vanishes and the statement is void. 
On the other hand, if /x > for all x G S 1 at t = 0, the right hand side 
of (29) is a positive number because S 1 is compact. The conclusion follows 
from the fact that if (r, s) G Us then 

p — V- — \{ r + s) — (p^ > tpo — (f- — 8 > 0. (30) 

□ 

4.2 The blowup equation 

For a homogeneous system of two conservation laws, it is possible to show 
that shock waves will form for suitable initial data by introducing a new 
unknown which is a rescaling of the spatial derivative of r. The scale factor 
can be chosen such that the new unknown satisfies a differential equation 
which does not involve derivatives of s. We follow a similar route, but as we 
will see below it is a bit harder to decouple the equations when source terms 
are present. 

We are interested in the spatial derivative ei(r). First note that 

[ei, e ] = c°i e + c^ei = c° 10 D + + (c\ - Kc° 10 )e 1 (31) 

from the definition of D + . It follows that 

D+ei = (ei - c° w )D + - (/c r ei(r) - re s ei(s) + c\o - Kc°i )ei. (32) 

We adopt the convention of denoting partial derivatives by subscripts, and 
we also write f ei = e\ f t + t\f x for the partial frame derivative of / with 
respect to e±, i.e., e\{f) with r and s are regarded as constants. 
Applying the operator D + ei to r and using (25a) gives 

D+ei(r) = -/t r ei(r) 2 - « a ei(s)ei(r) - - /tc°i )ei(r) 

+ / ei + /rei(r) + / s e 1 ( S )-c ( W. (33) 

In particular, there is a term quadratic in ei(r) which might cause blowup of 
ei(r). The problem is that there is a term involving both e\(r) and ei(s). As 
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in the homogeneous case, the mixed term can be eliminated by introducing 
a function h(r, s) by 

h(r,s)= / (K-A)- 1 -d,, (34) 

J 7 (r) OS 

where 7 is defined in Lemma 4.1. Now let R = —e h ei(r). Differentiating R 
and simplifying gives 

D+R = a 2 R 2 + ai-R + a - e h f s e 1 (s), (35) 

where 

a 2 = e~ h K r , (36a) 

Ol = /r + //i r + (k - A) _1 #k s + C01 - KC°oi, (36b) 

ao = -e /l (/c° 01 + / ei ). (36c) 

We would like to use (35) to show that R has to blow up for appropriately 
chosen initial data. But there is still a problematic term involving ei(s) which 
cannot be estimated using bounds on r and s alone. This term is not present 
in the homogeneous case since then f s vanish. It is of course possible to 
define a quantity similar to R based on ei(s) instead of e.\{r). The result is 
a system of two coupled equations, but then we cannot use ODE techniques 
directly. 

It is, however, possible to remove the ei(s) term by the following proce- 
dure. By construction, 

D + s = D_s + (k - A)ei(s) =# + («- A)ei(s), (37) 

where we have used (25b). Since k — A 7^ 0, we can introduce a function 

<j>(t,x,r,s) = - f e h (n — \)~ 1 f s ds. (38) 

J 7 (r) 

By the chain rule, (37) and (25a), 

D + (j) = eo + K(p ei + 4> s (k - A)ei(s) + <p s g + <f> r f, (39) 

where eo and ei are the frame partial derivatives of 0, regarding r and s 
as constants. Applying (38) to (39) gives 

D + <p = eo + K<p ei - e h f s ex (s) + <p s g + <p r f. (40) 
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Now (35) may be written 

D + (R - ( f ) ) = A 2 (R- 0) 2 + A 1 (R -(/>) + A , (41) 

where the coefficients are 

A 2 = e~ h K r , (42a) 

M = fr + fh r + (k- \y l gn s + Coi - KC°oi + 2e~ h K r (j), (42b) 

A = e- h n r 4> 2 + (/ r + //i r + (k - A)~ V* + cV - kc° O i)0 (42c) 

- e fc (/c° i + / ei -/.(«- A)- 1 ^) - r / - eo - K <f> ei . (42d) 

Note that (41) is a first order ODE in R — whose coefficients can be esti- 
mated in terms of r and s. 



4.3 Blowup of i? 

We start with a simple lemma about blow-up of solutions to an ODE with a 
quadratic nonlinearity (see, e.g., [1], p. 72 or [9], Lemma 1.3.2). 



Lemma 4.4. Let v(t) be a solution on [t ,ti) of 
dv 



dt 



A 2 (t)v 2 +A 1 (t)v + A (t), v(t )=v , 



(43) 



where A 2 , Ai and A are continuous and bounded on [t ,ti) with A 2 > 0. 
Put 



K 1 {t)= j A 1 dt, K (t)= j |A |e- Xl dt, K 2 (t)= j A 2 e K 'dt. 

J to J to J to 

If Vq > K (ti) then 



K 2 (h) < (vo-Koih)) 1 
and we have the following lower bound on v in [to,ti): 

e- K ^v(t) > K (h) - K Q (t) + \{v - Mt^y 1 - K 2 (t) 



(44) 
(45) 



-l 



(46) 



Proof. The linear term can be dealt with by putting V = ve Kl , which 
transforms (43) into 



^ = A 2 e K W 2 + A e~ K \ V(t ) = v , 
dt 



(47) 
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Let W be the solution of 

^ = A 2 e Kl (W - Koih) + K (t)) 2 - \A \e- K \ W(t ) = v , (48) 

which can be found explicitly as 

W(t) = K (h) - K (t) + [{v - Mt^y 1 - K 2 (tj\ \ (49) 

We want to show (46), i.e., V > W. Now 

^ < A 2 e K ^W 2 + A e~ K \ (50) 
at 

so subtracting (47) from (50) gives 

^-(W-V) < A 2 e Kl (W + V)(W -V), (51) 



and a Gronwall estimate implies that W — V < as long as W + V > 0. 
Whenever W - V < 0, W + V > 2W, so V > W as long as W > 0. 
By definition, W > when K 2 (t) < (v - ^(ti)) -1 - But W -> oo as 
-^(i) — > (fo — ^o(ii)) -1 , so since F is bounded on [io,*i) we must have 
K 2 (t 1 )<(v -K (t 1 ))- 1 . □ 

Note that in [1] and [9], the factors involving K 1 are moved outside the 
integrals. We avoid this since we want to keep the simple form of (46). 

Lemma 4.5. Let U be an open set in R 2 and put Q = [to,ti) x S 1 x U. 
Assume that 

1. (r, s) remains in U for t G [to,ti), 

2. e ° + nei° is continuous and positive on Q, 

3. A 2 , A\ and A , given by (42), are continuous with A 2 > and 

C2 = infA 2 , C\ — sup|Ai| and Co = sup|A | (52) 
u q n 

are all finite with C 2 > 0, and 

4- 4> is bounded on Q. 
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// 

R(t ,x) > <j)(t ,x,r (x),s (x)) + Coxih) + (C 2X (ti)) _1 + Ci/C 2 (53) 
for some x G S 1 , where 

< {e c^_ l)/Cl !/Cl>0 , 

!/C 1= 

and £(£i) is t/ie parameter value along the characteristic g x corresponding to 
the time t\, then R cannot be bounded on [to,ti] x S 1 . 

Proof. Assumption 1 ensures that the other assumptions apply along the 
characteristic g x for t G [io^i]- Along g x , the parameter £ satisfies dt/d£ = 
eo° + /cei°, so it follows from assumption 2 that £ is a continuously differen- 
tiable and increasing function of t along g x . 

Because of assumption 3 and assuming that R — (p is bounded, we can 
apply Lemma 4.4 to (41) along g x , giving 

R(to,x) -(j>(to,x,ro(x),s (x)) < tf (£(*i)) + ^(h))' 1 . (55) 
If C\ > 0, we have 

|^| < Ctf(f), < ^(e^-l) and X 2 > ^(1-e^)), 

(56) 

so 

i2(*o,x) - 0(t o ,x,r o (x), s (a:)) < C„x(*i) + (C^i))^ + C1/C2, (57) 

where x(*i) = (e Cl ^^^ — 1)/Ci. When Ci = the same inequality holds but 
with x(ti) = €(ti)- This contradicts the assumptions of the lemma and so 
in fact R — is unbounded. It follows that R must blow up since does 
not. □ 

Theorem 4.6. Let Us and fls be as in (26). Assume that 

I. A ; A r; \ s , K r and k s are continuous and bounded on Us, with n r and 
k — A bounded away from 0, 

II. /, f r , f s , f seo , f ei and g are continuous and bounded on 
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III. e ° + nei° and e ° + Aei° are continuous, bounded, positive and bounded 
away from on an d 

IV. c°oi and c\\ are continuous and bounded on [to,ti) x S 1 . 
Suppose also that 

ei(r)(t ,x) > e- /l ^' s »)(0(t o ,a;,r o , So )+Cox(t 1 ) + (C 2 x(ti)) _1 +C 1 /C 2 ) (58) 

for some x G S 1 , where ti = min{t + 2~ 1 / 2 5EF~ 1 , t{\, F and E are given 
by (28), h is given by (34), by (38), Cj by (52) and x by (54). T/ien £/jere 
no smooth solution of (25) on [t ,ti\ with initial data (r ,s ). 
Moreover, the right hand side of (58) can be estimated in terms of \\ro\\, 
\\ s o\\> ^0; ti, x , ?> an d the bounds in assumptions I-IV. 

Proof. First of all, that (r,s) G Us for t G [to? ^i) follows from Lemma 4.2, 
which applies because of assumptions II and III. Thus condition 1 of Lemma 4.5 
holds. Condition 2 follows directly from assumption III. 

To simplify the presentation we will denote the supremum norms over 
S 1 , Us and fl s with the same symbol ||-||. From the definition (34) and 
Lemma 4.1, h is continuous on Us with 

\\h\\ < ||(/€-A)- 1 ||||^||(2|| So || (59) 

It follows that h is bounded because of assumption I. Differentiating (34) by 
r and performing a partial integration to get rid of the mixed second partial 
derivative n rs gives an estimate 



HM < ||(k- A)- 1 ||(2||K r || + ||kJ) 

+ ||( K -A)- 1 || 2 (|| Kr ||||A s || + ||^||||A r ||)(2|| So ||+5), 



so it follows from assumption I that h r is bounded as well. 

From the definition (38) and assumptions I and II, <p is continuous on Vts 
with 

||0||<||(K-A)- 1 ||elHI||/ s ||(2|| So ||+5). (61) 

Differentiating with respect to r and performing a partial integration to 
avoid the term involving f rs gives an estimate 

110,11 < e l|fe " || (k - A) _1 || [2||/ r || + \\f s \\ + (HMIHifJ + \\f s \\\\H r \\)(2\\s \\+5) , 

(62) 
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where H = h — ln(/t — A). Differentiating with respect to e gives the 

6 S 1 1 1X1 db tj 6 

||0 eo ||<eM||( K -A)- 1 ||||/ se() ||(2|| So ||+5) ) (63) 
while for <f> ei we perform a partial integration before estimating, giving 

||0 ei ||<eW||( K -A)- 1 ||||/ ei ||(2+ 11^11(211^11+5)). (64) 

By assumptions I and II, 0, r , eo and <f> ei are bounded. We conclude that 
all terms in (42) are continuous and bounded. 

Now infu s A 2 > because of assumption I, so all conditions of Lemma 4.5 
are met. Since R = e h ei(r), ei(r) cannot be bounded on the characteristic 
g x if (58) holds. 

Finally, the only quantity in the right hand side of (58) whose dependence 
on the estimates has not been established by the arguments above is 
From the definition (28) of E and F, t\ can be estimated both from below 
and above by to, t\, x, 5, ||r ||, ||so|| an( i the bounds in assumptions I-IV. 
That the same holds for x(ti) follows from the definition (54) of x together 
with the estimate of A 1 above. □ 

The condition on f seo can of course be replaced by a condition on f eo by 
performing a partial integration as we did with h r , f rs and f sei . The reason 
for not doing so will become evident below. Note also that there is no explicit 
dependence of f eo in (42). 

Let us define the blowup time as 

= sup{t > t ] (r, s) is smooth on [to^]}- (65) 
Then we also have the following. 

Corollary 4.7. Under the conditions of Theorem 4-6, the blowup time t* 
satisfies x(t*) — X* where x* is the smaller root of 

ei (r)(t ,x) =e-^^(0(t o ,x,r o , So ) + C o x+ (C 2X ) _1 + C^) . (66) 

5 Blowup in spacetime 

We will now combine the results of the preceding sections to show finite life 
span of solutions to the Einstein-Euler equations. We first need to control 
higher derivatives of some geometric quantities, given bounds on the matter 
variables. 
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Lemma 5.1. If the conditions of Lemma 3.1 are satisfied and r and s are 
bounded with r + s > 2ip_, the following quantities are bounded: 



K',A",a",d t K,(3" and d t A' . 



(67) 



The bounds involve only those in (6) and the a priori bounds on r and s. 

Proof. From Lemma 3.1, p, j and S can be bounded in terms of r and 
s for a given equation of state. Using these bounds together with (5) in 
(4a), (4b), (4c) and (4f) immediately gives bounds on K', A", a" and dtK. 
Differentiating (4d) by x gives (5" expressed in a, a', K, K' and t, so it is 
bounded too. Finally, differentiating (4e) by x gives a bound on d t A' in terms 
of the previously obtained bounds. □ 

The previous results can now be combined to show a blow up result for 
the Einstein-Euler equations in plane symmetry. 

Theorem 5.2. Lett < ti < and suppose that the conditions of Lemma 3.1 
hold. Then there are smooth initial data on the surface t = t such that the 
corresponding smooth solution of the Einstein-Euler equations (4) and (11) 
does not extend beyond t — t±. 

Proof. Theorem 2.1 shows that all the geometric quantities in (5) are bounded 
in terms of the quantities in (6). Let Us and Qs be as in (26), with 5 as in (29). 
Then Us and fls depend on ro, so, to, t\ and 5, and 5 is limited by y?o — V 9 -) 
which in turn depends only on the equation of state and the minimum of 
r + s . By construction, r and s can be bounded on in terms of 5 and 
the maximum values of r and Sq. It then follows from Lemma 5.1 that the 
geometric quantities (67) also are bounded, in terms of (6) and 



We need to establish the bounds in Theorem 4.6. First Theorem 2.1 and 
Lemma 5.1 give that 



<*> ll r oll, IK 



(68) 



e ° + ne? = e ° + A ei ° 



= a 



-l 



(69) 



are continuous, positive and bounded away from and 



c° i = a 1 ei(a) and c\i = a 



1 Ae 1 ((3)-A- 1 e (A) 



(70) 



18 



are continuous and bounded in terms of (6). 

From the definition of u and the construction of Us, ln(l — u 2 ), (1 — u 2 )^ 1 
and (1 ± uw)" 1 can be bounded in terms of (68). Also, 

1 1 

- inf{r + s } - 5 < ip < - sup{r + s } + 5, (71) 
Z S 1 Z gi 



so for a given equation of state, Lemma 3.1 gives that w 1 , (1 — w 2 ) 1 and 

dui 

dip 



" : " are bounded in terms of (68) and 



(V? -<p-- S)- 1 = inf{r + s } - <P- ~ <*) (72) 

The balance law coefficients can be read off directly from (16). They are 

, , u + w 

« (r ' a) = TW (73a) 

A(r,a) = f^-, (73b) 
1 — «w 

tw + ^M-2MwA' 1 ei(A) _j 
fyt.x.r.s) — a ei a , 73c 

tw — Ku — 2uwA~ 1 eAA) , . . 

g(t,x,r,s) = ha ei(m. (73d) 

1 — Mtt) 

From the previously obtained bounds it follows that k and A are continuous 
on Us and can be bounded in terms of (68). Differentiating and using that 
w r = w s = |^ and u r = —u s = |(1 — w 2 ) we get 

2(1 + ww)^ V dy?/ 

with similar expressions for /c s , A r and A s . These quantities are continuous 
on Us and can be bounded in terms of (68) and (72). Moreover, n r and 

K _ A .Mi^) (75) 

1 — U 2 W 2 

are positive, and k" 1 and (« — A) -1 are bounded in terms of (68) and (72). 

From (73), / and g are continuous on Vis and bounded in terms of (6) 
and (68). Since w r = w s = |^ and u r = —u s — |(1 — w 2 ), the same holds 
for their r and s derivatives. 
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Finally the quantities appearing in f seo and f ei not present in / are d t K, 
K', d t A', A" and a". It follows from Lemma 5.1 that / seo and f ei are con- 
tinuous and bounded in terms of (6), (68) and (72). Note that f seo does not 
depend on d t a' since the term in (73c) containing a is independent of s. This 
is the reason for not replacing f seo by f eo in Theorem 4.6. 

Now all the conditions of Theorem 4.6 have been established and we can 
conclude that if the initial data satisfies (58) for some x, there is no smooth 
solution on [to, ti] . By definition, ei(r) = A~ 1 d x r. Since all quantities in (58) 
except d x r can be estimated in terms of the geometric initial data (6) and 
the matter initial data (68) and (72), choosing initial data with sufficiently 
large gradient d x r will prevent a smooth extension up to t — t\. 

□ 

The question, how big the class of solutions of the constraints is which 
satisfy the hypotheses of the theorem will not be treated in great generality 
here. It will, however, be shown that there are some data which do so for 
which the time of existence is arbitrarily short. Consider data for which 
u is identically zero. Then a particular solution of equation (4a) is given 
by K = t/3. Substituting this into the remaining constraint (4b) gives 
(VA)" = — |A 5 / 2 (— | + lQnfi). The equation to be solved is a special case 
of one treated in [10] where the existence of a solution A was shown by the 
method of sub- and supersolutions. This method also gives uniform bounds 
for the solution which in particular shows that the derivative // at a point 
can be made arbitrarily large within a family of solutions while maintaining 
a fixed L°° bound for A. At the same time fixed L°° bounds for /i and 
fi^ 1 can be maintained. This suffices to show that e\{r) becomes arbitrarily 
large within the family and provides the desired data for which the time of 
existence is arbitrarily short. 

6 An extension result 

In the previous section we established that if the initial data has sufficiently 
large gradients, there is no smooth solution of the Einstein-Euler equations 
beyond a certain time. It seems quite plausible that the obstruction to ex- 
tending the solution is the blowup of first derivatives of the matter variables, 
although it must be pointed out that we have not shown that. To inves- 
tigate this further, we will show that if the first derivatives of the matter 
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variables are bounded then the spacetime can be extended. The argument 
will be similar to that in [12]. First we need bounds on r and s, which can 
be established under a mild extra condition on the equation of state. 

Theorem 6.1. Assume that in addition to the conditions in Lemma 3.1, 
dp/dfi is bounded away from 1. Then the following quantities are bounded: 

r,s,p,j,S. (76) 

Proof. The extra condition and Lemma 3.1 implies that the factors in the 
denominators of (73) are positive and bounded away from 0. Thus the right 
hand sides of (25) can be bounded in terms of the geometric quantities (5). 
Integrating (25) shows that r and s are bounded along the characteristics. 
We may pass from these estimates to estimates in time because of (69). 
Finally, p, j and S must also be bounded because of Lemma 3.1. □ 

Lemma 6.2. Under the conditions of Lemma 3.1, if all spatial derivatives of 
order up to n of the geometric quantities (5) and the matter quantities (76) 
are bounded, the spatial derivatives of order up to n + 1 of the geometric 
quantities (5) are bounded. 

Proof. Note first that some of the spatial derivatives are already included 
in (5). As was shown in the proof of Lemma 5.1, we can solve (4) for the 
remaining quantities a", A", d t A', K' and (3". Differentiating the resulting 
equations n times gives the desired bounds. □ 

Lemma 6.3. Suppose that the conditions of Lemma 6.2 hold, dp/d/i is 
bounded away from 1 and d x r and d x s are bounded. Then the spatial deriva- 
tives of the matter quantities (76) of order up to n + 1 are bounded. 

Proof. Note first that because of the definition (2) of ei, using e" instead of 
introduces spatial derivatives of A^ 1 of order at most n — 1, which are 

bounded by assumption. 

Put r n — e"(r) and s n — e"(s). Applying (32) to r n+ \ gives 

D+r n+1 = (ei - c° 10 )D + r n - (/c r ei(r) - re s ei(s) + c\ - Kc° w )r n+1 . (77) 

If we can show that (77) is a linear equation in D + r n+ i with the coefficients 
given by smooth functions of spatial derivatives of r and s of order up to n 
and of the geometric quantities (5) of order up to n + 1, we would be done 
because of Lemma 6.2. 
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Putting n = 1 and using the expression (33) for D + r\ shows that D + r 2 is 
linear in r 2 , and the coefficients have the right dependence because of (73) 
and (70). For n > 2 the coefficients of D + r n are differentiated by x at most 
once in (77) so the result follows by induction. The case with s is completely 
analogous. □ 

Lemma 6.4. Under the conditions of Lemma 6.3, if all derivatives of the 
geometric quantities (5) and the matter quantities (76) of the form d\d™ with 
n > and < k < m are bounded then the <9™ +1 <9™ derivatives of the matter 
quantities (76) are bounded. 

Proof. As was the case with e\ and d x , we may use e instead of d t because 
of (2). The balance laws (25) can be used to bound an extra time derivative 
by introducing an extra space derivative, which is bounded by Lemma 6.3. 

□ 

Lemma 6.5. Under the conditions of Lemma 6.3, if 

1. all derivatives of the geometric quantities (5) and the matter quanti- 
ties (76) of the form <9 t fc <9™ with n > and < k < m are bounded, 
and 

2. all derivatives of the matter quantities (76) of the form d™ +1 d™ are 
bounded, 

then the derivatives of the geometric quantities (5) of the form <9™ +1 (9™ are 
bounded. 

Proof. See [12]. □ 

Putting together Theorem 6.1 and Lemma 6.2- 6.5 we see that under 
the hypotheses of Lemma 6.3, all derivatives of the geometric and matter 
quantities are bounded. When all the derivatives are bounded on a given 
interval, the solution can be smoothly extended to the closure of that interval. 
This results in a new initial data set and applying the local existence and 
uniqueness result from [12] gives an extension of the solution. We thus have 
the following theorem. 

Theorem 6.6. Let a solution of the Einstein equations with plane symmetry 
be given. Suppose that when coordinates are chosen which cast the metric into 
the form (1) with constant mean curvature time slices, the time coordinate 
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takes all values in the finite interval [to,ti) with ti < 0. Assume that the 
equation of state satisfies the conditions in Lemma 3.1, and that dp/dfi is 
bounded away from 1. If the first spatial derivatives of the energy density 
\i and the velocity parameter u are bounded on [to,ti), the solution can be 
extended beyond t — t±. 

7 Discussion 

In this paper it has been shown that plane symmetric classical solutions of 
the Einstein-Euler system can break down in finite time under some general 
assumptions on the equation of state. Time here is measured with respect to 
a CMC foliation. That the formation of singularities is due to the fluid and 
not just to a breakdown of the CMC foliation is shown by the following fact. 
If instead of the Euler equations we describe the matter content of spacetime 
by a collisionless gas satisfying the Vlasov equations and maintain all other 
assumptions the solutions exist globally in the future [12]. In the case of the 
Euler equations for a fluid without pressure (dust) a non-existence result was 
proved in [13]. The proofs of all these results have common elements. It is 
shown that the spacetime geometry has certain good properties right up to 
the singularity. Using this it is shown that when there are no singularities 
in flat spacetime there are also none after coupling to the Einstein equations 
(Vlasov) while in the cases there are singularities (Euler with pressure, dust) 
the mechanisms of blow-up in the flat space case (formation of shocks, blow- 
up of the density) provide a reliable guide as to what quantities are bounded 
or not and how to prove it in the case of coupling to the Einstein equations. 

A family of data along which the time of existence becomes arbitrar- 
ily small can be used to show that the results of [10] on solutions of the 
Einstein-dust system generalize to the case of fluids with non-zero pressure. 
The conclusion is that there are data for the Einstein-Euler equations with 
an equation of state of the kind considered in the present paper on a CMC 
hypersurface such that the corresponding maximal Cauchy development can- 
not be foliated by compact CMC hypersurfaces. To see that the proofs of [10] 
can be adapted note first that under the assumption u — the constraints 
for the fluid with pressure are identical to those for dust so that many of the 
arguments can be taken over directly. Cauchy stability for the Einstein-Euler 
system with < * < 1, which is necessary for the argument, follows from 
the fact that it can be written in symmetric hyperbolic form so that standard 
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results apply. See for instance [16], in particular p. 140. The final element 
which is required is to know that the argument for breakdown of the solution 
can be localized in space. Since the contradiction to existence on a certain 
time interval is obtained by an analysis of the evolution of certain quanti- 
ties along a single characteristic it is only necessary to show that in a short 
time interval the characteristic can only move a short (coordinate) distance 
in space. This follows immediately from the available a priori bounds for the 
geometry. 

The assumptions on the equation of state in Theorem 5.2 are compati- 
ble with cases of physical relevance. Consider, for instance, the relativistic 
polytrope with 

fi = m + Knm 1+ «, p = Km 1+ « (78) 

where n > 1 and K > are constants and m the rest mass density. In this 
case all the assumptions are satisfied. They are also satisfied for an equation 
of state leading to a linear relation p = Kfi with K > 0, as considered 
in [17] and [3]. This is true in spite of the fact that if a full thermodynamic 
treatment of this case is attempted the internal energy is negative at low 
densities and so its physical interpretation is questionable. 

In the case of Theorem 5.2 the density remains bounded on the maxi- 
mal interval of existence of a classical solution. This cannot be expected to 
hold for dust, although it is not actually proved in [13] that the density is un- 
bounded. Thus some restrictions on the equation of state are to be expected. 
It is probable that blow-up of the density occurs when the equation of state 
is such that the pressure is a bounded function of the density (cf. [19]). In 
Theorem 5.2 this case is excluded by condition 4 of Lemma 3.1. 

Ideally, there would be a useful interaction between the results of [3] 
and those of the present paper. It would then be possible to conclude that 
the classical solutions with finite-time breakdown can be continued as weak 
solutions and on the other hand that the weak solutions which have been 
constructed are not always classical. Unfortunately the time coordinates used 
in the two cases are different so that a direct comparison is not possible. 

Finally, some possible extensions of the results of this paper will be men- 
tioned. In the case of a special relativistic fluid the analogue of Theorem 5.2 
holds with the same proofs. In that case the time is the ordinary Minkowski 
time. Also, Theorem 6.6 shows that under the extra condition that dp/dfi 
is bounded away from 1, a solution for which the first derivatives of m and 
u remain bounded remains smooth, i.e. it can be extended smoothly to a 
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longer time interval. The results obtained here therefore support the picture 
that the breakdown of the classical solutions is due to the formation of shock 
waves. 
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